Phenomenological study on the significance of the scalar potential and Lamb shift 
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We indicated in our previous work that for QED the contributions of the scalar potential which appears at the 
loop level is much smaller than that of the vector potential and in fact negligible. But the situation may be 
different for QCD, one reason is that the loop effects are more significant because a s is much larger than a, 
and secondly the non-perturbative QCD effects may induce the scalar potential. In this work, we phenomeno- 
logically study the contribution of the scalar potential to the spectra of charmonia. Taking into account both 
vector and scalar potentials, by fitting the well measured charmonia spectra, we re-fix the relevant parameters 
and test them by calculating other states of the charmonia family. We also consider the role of the Lamb shift 
and present the numerical results with and without involving the Lamb shift. 

PACS numbers: 1 1.10.St,12.38.Aw,12.39.Pn 



I. INTRODUCTION 

By a symmetry consideration, Chen et al. suggested fH] that 
for the Coulomb interaction, to maintain the hidden symmetry 
SO(4) in the Schrodingger equation, the scalar and vector po- 
tentials must have the same weight in the Dirac equation. The 
hidden symmetry is just the familiar Lenz symmetry which 
also exists in the classical physics. However, if so, the orbit- 
spin coupling would disappear. In fact the scalar and vector 
potentials make opposite contributions to the orbit-spin cou- 
pling, thus if they have the same weight, their contributions 
would exactly cancel each other. It definitely contradicts to 
the data. Therefore, one concludes that this symmetry does 
not exist in the relativistic extension. Usually, one is tempted 
to think that the relativistic Dirac equation should possess a 
higher symmetry than its non-relativistic approximation, but 
this is not the case we are confronting. A general theory 
which only considers the Lorentz structure of the vertices [2], 
there are five types of coupling. But which one dominates 
should be selected by the underlying physics. We turn to look 
at the deeper side, namely start to investigate the problem in 
the quantum field theory. 

The basic theory which induces the electric Coulomb po- 
tential is QED whose coupling is vector-type ^7 yi\i , thus 
at the tree-level, the induced potential is the vector one and 
the other types should be induced at higher order, i.e. loop 
level. In our earlier work [3], we showed explicitly that the 
scalar coupling 101 which results in the scalar potential, ap- 
pears at the loop level and its contribution is suppressed by 
a factor a/ir. For QED it is a small value and cannot make 
a sizable contribution. Thus the apparent SO(4) symmetry at 
the classical level is almost fully violated. However, the sit- 
uation would be different for the QCD case, because first a s 
at the charm-mass-scale is much larger than a and secondly 
the non-perturbative QCD effects may also cause the scalar 
potential. 
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This case is noticed by Leviatan and some studies have been 
carried out [4-6]. In this work, we are not going to further dis- 
cuss the origin of the scalar potential or try to derive it from 
the quantum field theory, but generally assuming its existence 
and by fitting the spectra of the charmonia family, we obtain 
its fraction. Moreover, the QED theory predicts the Lamb shift 
which is due to the vacuum effects. In QM, it only shifts the 
S-wave spectra because in the non-relativistic limit, it is pro- 
portional to <5(r), but by the quantum field theory, the other 
Z-states are also affected. In other words, by considering the 
Lamb shift, the positions of the spectra would deviate from 
that obtained without the Lamb shift. In this work, we in- 
clude its contribution and re-fit the charmonia spectra to ob- 
tain a new set of the model parameters. For a comparison, we 
will present the numerical results with and without taking the 
Lamb shift into account. 

Unlike the hydrogen-like atoms where the nucleus is very 
heavy and approximated at rest, therefore only the motion of 
electron is considered and the corresponding equation, either 
the relativistic Dirac equation or non-relativistic Schrodinger 
equation, is a one-body equation. However, for charmonia, 
the charm and anti-charm quarks are of the same mass and 
the equation which properly describes charmonia, should be a 
two-body equation. 

For simplicity but without losing the significant characters, 
we do not directly solve the two-body Dirac equation which 
is very complicated. One can derive the effective potential 
between the two constituents (c and c) in terms of the pertur- 
bative theory where the effective Lorentz vertices are set ac- 
cording to the general Lorentz structures [7]. Because of the 
limitation of the perturbative theory, we can only obtain the 
Coulomb-type interaction and the corresponding spin-orbit, 
spin-tensor and relativistic correction pieces. It is noted that 
the fundamental QCD indeed provides only the vector poten- 
tial at the tree level, but as indicated above, the loop effect 
and even non-perturbative effect may result in scalar poten- 
tial. Thus we just keep the potential forms and introduce two 
phenomenological constants in front of the scalar and vector 
potentials and the induced terms. For the confinement piece, 
we employ the linear confinement i.e. the Cornell-type. In 
fact, the exact form of the full potential including both scalar 
and vector pieces was given by Lucha et al. |7[] and we just 
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re-check their results and then substitute the potential into our 
Schrodinger equation. 

Now we can reduce the two-body Schrodinger equation into 
one particle equation where the kinetic term is j^p 2 where fi 
is the reduced mass and is m c /2 in our case. Solving the 
differential equation, we obtain the spectra. Since there exist 
several phenomenological parameters which so far cannot be 
derived from the underlying theory, we can fix them by fitting 
a few well measured charmonia states. 

Moreover, as well known, the vacuum fluctuation induces 
the Lamb shift. The basic Lagrangian of the Lamb shift has 
been derived by some authors, and for interaction, we have 
Hint = —Lint II HI]- Thus we substitute the expression 
into our data fitting process to re-derive the phenomenologi- 
cal parameters. Indeed, the Lamb shift only occurs at the loop 
level, but the Coulomb-type —a s /r appears at the tree level 
of QCD. It seems that they belong to different levels, but as 
we introduce the phenomenological parameters which include 
the loop and non-perturbative QCD effects, we cannot distin- 
guish between the tree level contribution and that of higher 
orders . However, for the Lamb shift, we do not introduce a 
new phenomenological parameter, but use the derived form 1 . 

There are some subtleties in the calculations which we will 
address in the text. 

This paper is organized as follows. In Section HI1 and HIT1 
we introduce the generalized Breit-Fermi Hamiltonian and the 
Schordinger equation for the cc bound states: J/ip, x c o(lP), 
Xci(lP), ?/c(2S) and ip(2S). Then we numerically solve the 
eigen-equations for these bound states and fix the parameters. 
In Section |IV] the Lamb shift is concerned and another set of 
the parameters is given to improve our predictions. The last 
section is devoted to our conclusion and discussion. 



II. THE GENERALIZED BREIT-FERMI HAMILTONIAN 
AND SCHORDINGER EQUATION 

For the cc meson, the generalized Breit-Fermi Hamiltonian 
was given in Refs. [2J, LlOQ as 



H = H + H 1 + ..., 

H Q = V — + S{r) + V(r), 
to 

Hi = H s( j + H si , 

H s d = His + H ss + H t 
1 
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-V'(r) + V" - rV 



(le) 



where, V and S stand for the vector and scalar potentials 
and H S i and H s d represent the spin-independent and spin- 
dependent pieces respectively. For the linear confinement 
piece we adopt the Cornell potential lfTTll . Thus the total po- 
tential at lowest order reads 



U(r) = V(r) + S(r) = -aC F — + b K 2 r 

r 



where, 



V(r) = -c CfoisIt + dn 2 r 

S{r) = -(a- c)C F a s /r + (b - d)n 2 



(2a) 



(2b) 



With the Hamiltonian ([]]) and the potential (f2]), one can solve 
the Schordinger equation 



{E - 2ro)*(r) = HV(r) = (H + ffi)*(r). 



(3) 



(4a) 



If we define the radial wave function as R(x) with the di- 
mensionless variable: x — nr, then the radial equation is writ- 
ten as 2 

where, 

A(x) 

with 



—ru(x) = A{x)u(x) 

dx z 



-m (E- 2m- U(x) - H x 
-^(E -2m-U(x)] 2 



1(1 + 1) 



(4b) 



to = m/re, E — E/k, 

Hi = Hi/ k, U(x) = U(x)/k. 



The approximation 



is used in ([3). 



m(E-2m- U(r)) 



(4c) 



(4d) 



III. THE ENERGY GAP FUNCTION OF THE cc 
CHARMONIA AND THE NUMERICAL RESULTS 

The radial equation dU can be solved in terms of the method 
called " the iterative numerical process " which is introduced 
in literatures, (for example, see lfl2j[T3ll ). We have improved 



It is noted that for a formula which is derived in the field theory, one can 
separate the contributions corresponding to different orders as long as there 
is no phenomenological parameters involved, and that is the case we deal 
with the Lamb shift, please see the text for details 



2 The standard form of the radial equation can be easily found in fl2ll . and 
the method to make it dimensionless is given in fTlll . 
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this method, and then fix the parameters a, b, c, d by fitting 
the well measured spectra of cc charmonia: J/ip, Xco(lP), 
Xci(lP). ry c (2S) and ^(2S). Instead of directly fitting the 
masses, we construct a series of relations which should be 
fitted: 



' m M2S)] - m [xci(lP)] = £ [2 3 Si] - £ [l^ 
m V(2S) - m [J/^(1S)] = £ [2 3 Si] - £ [l 3 Si 
fo(2S)] = £ [2 3 Si] - £ [2%] ; 
[Xco(lP)] = E [2 3 Si] - £ [1 3 P ] 



to [^(2S)j 



m 
m 



(5) 



where, _E |n^ s+1 lj] represents the eigen-values of the radial 
equations @ with various quantum numbers n r , j, 1, and s. 
Because the parameters a, b, c and d are involved in the po- 
tential (f2j), E [n^ s+1 lj] must be functions of these parameters. 
to [meson] are the masses of the individual states which are 
shown in the following table lfl4ll : 



TABLE I: The experimental central values of the spectra of the cc 
charmonia states 



state 


m(GeV) 


state 


m(GeV) 


state 


m(GeV) 


JMi a S!) 

Xco(l 3 Po) 


3.0969 
3.4148 


Xd(l a Pi) 
»7c(2 1 So) 


3.5107 
3.6370 


^(2 3 Si) 


3.6861 



Sequentially, the parameters a, b, c and d are obtained by 
solving Eqs.(f5]l. By means of the Newton's iterative method, 
we have achieved as (Th e details about the numerical method 
can be found in Ref. [15].) 

a = 1.1715, b = 1.2250, c = 0.8087, d = 0.5291 (6) 

Here we set a s ~ 0.36 and k = 0.42 GeV which seem 
somehow different from the values given in literature fl6l - [l8ll . 
But as noticed, the deviation may be included in the phe- 
nomenological parameters a, b, c and d. 

A few words are about our choice of the input. In prin- 
ciple, any five well measured states of the charmonia can 
be used as the input. However, unfortunately, the relation- 
ship between the E [n^+Mj] and the parameters in (f5]l is 
complicated, taking the central values of the masses of J/ip, 
Xco(lP)> Xci(lP)) ?7c(2S) and ip(2S) as the inputs one can 
obtain reasonable solutions, otherwise, the equations (|5j do 
render solutions for a, b, c and d. The reason is due to the 
experimental errors. 

Given a, b, c and d in ©, the masses of the charmonia states 
can be written as: 



M 1 (n? + H j )=E[n? + H i ]+E Q 
where, Eq is the zero-point energy: 

E = m[J/i/)] - £[l 3 Si] 
and the final result is shown in Table-HIlbelow: 



(7) 
(8) 



TABLE II: The mass spectra for the charmonia states (in GeV), 
with m c = 1.84GeV. The mass of the EXP is the value given in 
PDGlil. 



meson EXP Prediction 


meson EXP Prediction 


77c(:PSo) 2.9803 3.0189 
J/^(l 3 S x ) at 3.0969 3.0969 
Xco(l 3 Po) flt 3.4148 3.4148 
Xd(l 3 Pi) fit 3.5107 3.5107 
M^Pi) 3 - 5259 3.5100 


Xc2(l il P2) 3.5562 3.5564 
J7 c (2 1 S ) flt 3.6370 3.6370 
</>(2 3 Si) flt 3.6861 3.6861 
V(3 3 Si) / 4.1164 



Explicitly, in the process, the masses of J/ip, Xco, "0(2 3 Si) 
and Xci(l 3 Pi) are taken as inputs to obtain the parameters 
and then the masses of other states in the family: ry c (lS), 
h c (lP), x C 2(lP) and -0(3S). The numbers are predicted. 



IV. THE MASS SPECTRUM AS THE LAMB SHIFT IS 
TAKEN INTO ACCOUNT 



As well known, the Lamb shift is due to the vacuum fluc- 
tuation and may cause sizable effects on the meson spectra. 
Indeed, the QED Lamb shift may not be very significant be- 
cause of smallness of the fine structure constant a, but for the 
QCD case, the situation will be different. 

On the other hand, the Breit-Fermi Hamiltonian used in 
Section-|In] does not include the effects of the Lamb shift in 
the eigen-energy (0]). In this section, we will take the Lamb 
shift into account. However, we do not introduce the Hamil- 
tonian induced by the Lamb shift into the differential equation 
because the corresponding pieces are very complicated and it 
is not necessary to do so. Instead, we simply add the estimated 
values of the effects to the binding energies of various states. 
Repeating the procedure done in last sections and adding the 
Lamb shift effects to the spectra, we re-fit the data to obtain a 
b c and d again and predict the mass spectra of the rest reso- 
nances. 

Namely, we set the mass of a bound state as: 



2m r 



AE 



LS 



M 



EXP 



(9) 



where E is the solution of the eigen-equation, AE-^g is the 
energy caused by the Lamb shift. Solving the equation, one 
can obtain the parameters again. 

The authors of Ref.iH, gave the theoretical expressions 
for the binding energies which involve contributions of the 
Lamb shift. When we only concern the Lamb shift, we must 
single it out from the general formulas. It is not difficult, as 
a matter of fact, because the Lamb Shift starts at 0(a 3 )|[ir 
The Lamb Shift can be written as: 



A£[n,j,l,s] 



m [A£(a 3 ) + 
+ AE(a 6 s ) + 



AE(a 4 s ) + AE(a 5 s ) 



(10a) 



For readers' convenience, let us directly copy Titard's for- 
mulas [9] below, where we dropped the tree-level terms and 
the relativistic corrections, and we have: 



4 



AE(a 3 s ) 
AE(at) 



C 2 

-a3^(2/3 o7£ +4 ai ); 

4 Cf // i A)s2 



7^(Mo + |) + (^+7l)f 



(10b) 
(10c) 



r 



Hoang et al. estimated the contribution of higher orders ically, these high-order terms can be attributed to the effects 
0{a 5 s ) and 0(a®) to the binding energies|8]. Phenomenolog- of the Lamb Shift: 



J 



AE(a 5 s 



AE(a e s 



5l Cf \C a 



C\ , AC A C F , 2C 2 



8 



2 n(21+l) n V 21 + 1 
4nl(l + l)(21+l) (4Xljs+ < 512 >ljs) / 5 



3n 

r 2 

a" log a 



47r 2 n 2 |_ 6n 
C A C F (1 - 6,0)6, 
8nl(l + 1 



36\pC F 
2n 



(CU + 2Cf) 



A)(^ - C F ) + ^(25C A + 22C F ) 



CaC f 6\q6 s i 



6n 



WTT) [l3o{2Xlis 



- < S 12 >lj S 



-C^(2Xij s + < Si2 >ij S )] 



(lOd) 



[5/3 + 7(7^] 



(10e) 



r 



n in Ref.dTOli stands for the principal quantum number as n = (also see B^ l2ll - l24tl ). 

n r + 1, where, n r and 1 are defined in Section-HITl All the The Lamb Shift A Mn,j, 1, s] depends on the coupling con- 
constants as ai, a,2, b\, Pi (i — 1, 2, 3) are given in Ref.[20] stant a s (see Eq.dTOll) [9] as: 



a s (p 2 ) 



2tt 



In A*/ A 



Aln(ln/x 2 /A 2 ) , fttfQnf/A^-ftlnQnif/A^-ft+fofo 



f3 2 In /i 2 / A 2 



ffiln 2 At 2 M 2 



(11) 



It is noted that unlike the others in the full Hamiltonian which 
can be written in the pure operator form, the contributions of 
the Lamb shift to the spectrum energies are always associated 
with the concrete states. 

Using the formulas given above, one can evaluate the Lamb 
Shift of the charmonia states. The scheme of renormalization 
is suggested by Pineda et al.ffl l2(il . Actually, there is a term 

2 



In 



in the theoretical expression of the energy (see 

where, a (a* 2 ) stands for the Bohr radius and fi is the 
renormalization scale: 

2 



«(At 2 ) 



where, 

Ma* 2 ) 



«i 



mC F a s (fi 2 ) 



7fA)\ ot 



(12a) 



IE 



A -l- & 
a-iPo + -5- 



(12b) 



If one defines [20] : 



A* = — 
na 



(13) 



this choice of At will cancel the terms related to ln[2|^] in the 
spectrum energy. 

The value of the parameter A is near 0.30 GeV[20]. Here, 
we choose it as 0.275 GeV. The reason is that, at this point, 
a" =2 = 0.38, just near to the value of a s we used in Section- 

m 

It is obviously different from the conventional renormaliza- 
tion scheme we commonly used. A consequence is that the 
coupling constant a s is different for different quantum num- 
ber n: 



0.31 , a r s 



0.38 , a 



n=3 



0.43 



(14) 



Simply adding the Lamb shift to the total binding energy 
is like that we change the zero-point energy for each state. 
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We still select masses of J/ip, Xco(lP)> Xci(lP)> ?7c(2S) and 
ip(2S) as inputs, and solve the equation (0 again as what we 
did in the last Section. But the value of a s in © is taken 
as that value given in Eq.(TT4b which depends on n. The new 
solutions of a, b, c, and d are: 



a W = 1.3943, = 1.4057, 
c« = 0.6243, = 0.9910; 

a( 2) = 1.4191, 6 (2) = 1.3292, 
c( 2 ) = 0.6459, d^ = 0.9438. 



(15) 
(16) 



where, expression ( TT3T > is the solution when the Lamb Shift is 
taken up to order 0(oP s ): 

AE[n,j,l, S }=m[AE(a 3 s )] 

and ( TToT l is for the Lamb Shift: 

AE[n, j,l, s] = m [A£(^) + A£(af) + AE(a 5 s ) + AE(a 6 s )] 

up to the order 0(af). With these two solutions, our predic- 
tions are given in Table-dlllb). 



TABLE III: The mass spectrum with the Lamb Shift (in GeV), where, 
the LS stands for the contribution of the Lamb Shift, Mi ,2 is the 
predicted mass when the parameter is set as in Eq,l|15t or Eq.l|16t 
and Mi, 2 stands for M{, 2 = Mi, 2 + A£ (1 ' 2) . 



meson 
^(LSo) 
J/^(l 3 Si) 
Xco(l 3 Po) 

Xcl(l 3 Pl) 
Xc2(l 3 P 2 ) 

?7c(2 1 S ) 
^(2 3 Si) 
^(3 3 Si) 



AE 



— 



M 2 



Mo AE 



— 



M 3 



Mi 



-0.0674 
-0.0674 
-0.0310 
-0.0310 
-0.0310 
-0.0310 
-0.0310 
-0.0310 
-0.020 



3.0820 
3.1643 
3.4458 
3.5417 
3.5593 
3.6079 
3.6680 
3.7171 
4.1460 



3.0146 
3.0969 
3.4148 
3.5107 
3.5283 
3.5769 
3.6370 
3.6861 
4.126 



-0.1196 
-0.1245 
-0.0799 
-0.0802 
-0.0803 
-0.0804 
-0.0714 
-0.0728 
-0.0531 



3.1612 
3.2215 
3.4946 
3.5909 
3.6063 
3.6552 
3.7084 
3.7589 
4.1746 



3.0416 
3.0969 
3.4148 
3.5107 
3.5260 
3.5748 
3.637 
3.6861 
4.1215 



the rest resonances of charmonia. It is shown that unlike the 
QED case where the fraction of scalar potential is very small 
and negligible, the fraction of scalar potential is of the same 
order of magnitude as the vector potential. This is consistent 
with the conclusion of Re f. 12511 and this is not surprising. As 
we indicated that for the vector-like coupling theories QED 
and QCD, the scalar potential can only appear at loop level or 
is induced by non-perturbative effect (QCD only). Thus it is 
loop-suppressed. However, for QCD, the coupling is sizable 
and the non-perturbative effects somehow are significant, so 
one can expect the fraction of scalar potential is large. 

Moreover, the Lamb shift is induced by the vacuum fluctu- 
ation and only appear at loop level, indeed the leading con- 
tribution is at 0(al). Therefore for the QED case, it is hard 
to observe the Lamb shift (observation of the Lamb shift is 
a great success for theory and experiment indeed), however, 
for QCD the effects are not ignorable. By taking into account 
the Lamb shift, we re-fit the model parameters and find they 
are obviously distinct from that without considering the Lamb 
shift. 

In the work, by studying the charmonia spectra we in- 
vestigate the contribution of higher orders of a s and non- 
perturbative QCD effects. However to distinguish between 
them, one needs to do more theoretical researches. This re- 
sult helps us to get a better understanding of QCD, espe- 
cially the non-perturbative effects. Even though it is only half- 
quantitative, it is an insight to the whole picture. 

When we take into account the contribution of the Lamb 
shift to the mass spectra, it is more obvious that higher order 
effects are important in QCD. Because the Lamb shift only 
appears at order 0(a^), its existence manifests higher order 
effects. Our calculations show that while higher orders up to 
0(a & s ) are involved, the fitted values of a,b,c,d are different 
from those when only O(a^) is considered. 

The same strategy can be applied to the bottomonia family 
and even the B c resonances where one can further test the 
theoretical framework and investigate the higher order QCD 
behaviors. That would be the contents of our next work. 



V. CONCLUSION AND DISCUSSION 

In this work, we study the role of scalar potential to the 
spectra of charmonia. Our strategy is that the scalar and vec- 
tor potentials have different fractions which manifest in their 
coefficients (in the text, they are a, b, c and d for the Coulomb 
and confinement pieces respectively). By fitting some mem- 
bers of the charmonia family, we can fit them. Then with the 
obtained parameters, we further predict the mass spectra of 
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